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We carry out a comprehensive study of the dynamics of large-scale perturbations in quintessence 
scenarios. We model the contents of the Universe by a perfect fluid with equation of state wt and 
a scalar field Q with potential V(Q). We are able to reduce the perturbation equations to a system 
of four first-order equations. During each of the five main regimes of quintessence field behaviour, 
these equations have constant coefficients, enabling analytic solution of the perturbation evolution 
by eigenvector decomposition. We determine these solutions and discuss their main properties. 
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I. INTRODUCTION 

Recent observations seem to indicate that the Uni- 
verse is undergoing a period of accelerated expansion [1] . 
Whereas cosmologists initially introduced a cosmologi- 
cal constant in order to explain this, a range of differ- 
ent models have emerged, amongst which quintessence 
has been particularly prominent in the literature [2, 3]. 
It is defined as a scalar field rolling down its potential 
and presently dominating the dynamics of the Universe. 
An important class of quintessence models are known as 
tracking models [2, 3], where the late-time evolution of 
the field has an attractor behaviour rendering its evo- 
lution fairly independent of initial conditions. In con- 
trast to a cosmological constant, which is by definition 
perfectly homogeneous, the quintessence field can, and 
indeed must, have perturbations. 

The evolution of perturbations in quintessence models 
have been studied by many authors [2, 4, 5, 6]. In this 
paper we carry out an exhaustive and elegant analysis 
of those in the large-scale approximation. We model the 
contents of the Universe by a perfect fluid with equation 
of state Wf and a scalar field Q with potential V(Q). We 
assume a flat Universe throughout. 



II. BACKGROUND EVOLUTION 

Before studying the perturbations, we recall some re- 
sults for the homogeneous background evolution. The ge- 
ometry of the Universe is described by a flat Robertson- 
Walker metric 



ds 2 = -dt 2 + a 2 (t)dx 2 . 
The Einstein equations 
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relate the matter components to the geometry. The in- 
dices "f" and "Q" always refer to the perfect fluid and the 



quintessence field respectively, and dots are time deriva- 
tives. We will use a prime to denote a derivative with 
respect to N = log(a/a ). 

The evolution of the fluid is straightforward, with its 
energy density scaling as a~ 3 ^ 1+W! \ where Wf is the 
ratio of pressure to energy density of the fluid. The 
quintessence field follows the Euler-Lagrange equation 



Q = -3HQ 



and its equation of state is 
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Depending on the precise model and on the choice of 
initial conditions, the quintessence dynamics can feature 
up to five main regimes, which were classified by Brax et 
al. [5] and which appear in a sequential order. During 
the first three the quintessence field is sub-dominant. 1 
The "kinetic" regime is characterized by the domina- 
tion of the kinetic energy which scales as a -6 . In the 
"transition" and "potential" regimes the potential energy 
dominates and the energy density remains constant. The 
sound speed of the quintessence field is defined by 
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and is equal to 1 or — 2 — Wf respectively during those 
regimes. During the "tracker" regime the quintessence 
field approximately mimics the behaviour of the fluid, 
and usually its energy density is still sub-dominant when 
tracking begins. If wq — wt there is perfect tracking, 
even if the quintessence field is not sub-dominant. Finally 
the field enters its "domination" phase, during which wq 
tends to —1 in most cases. 



1 Actually, one can consider initial conditions with domination of 
the quintessence field, but they usually lead to an "overshoot" 
of the required present energy density, and are therefore not in- 
teresting. 
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FIG. 1: An example of the evolution of the four parameters 
Qq, wq, CgQ, and 6q/5 in a realistic Universe, showing the 
five different regimes and also the transition between radiation 
domination and matter domination (at N ~ —7). We used 
an inverse power-law potential V(Q) = Vo(Q/rrip\)~ a , with 
a = 1 and Vb = 3 x 10~ 124 m%\. In this case we have the usual 
sub-dominant tracker regime. 



Four useful parameters can be defined to describe these 
five regimes. They are the quintessence density parame- 
ter Qq, the equation of state conveniently parametrized 
as 7q = 1 + wq, the speed of sound c 2 q as defined in 
Eq. (6), and one further parameter relating to the speed 
of sound defined as 
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In order to simplify the notation, we define the vector 
x = (ftQ,-fQ,c^Q,0Q). In Fig. 1, we show an example 
of the evolution of those parameters. We take the re- 
alistic case Wf = 1/3 during radiation domination and 
Wf = during matter domination, and we use an inverse 
power-law quintessence potential. We clearly see the five 
different regimes, and also the transition between radi- 
ation domination and matter domination which in this 
case occurs during the tracking regime. In Table I we 
give the values of the parameters in the general case for 
each regime. 
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TABLE I: Values of the four parameters Q.q, jq, c^q, and i 
during the different regimes. 



III. PERTURBATION EVOLUTION 

In order to describe the perturbations we choose the 
conformal Newtonian gauge [7]. As long as there is no 
anisotropic stress, the perturbed metric is 



ds 2 
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where in this case $ is equal to the gauge-invariant po- 
tential defined as in Ref. [7]. We work in Fourier space 
and compute the first-order perturbed Einstein equations 
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$ + 4ff$ + (2H + 3i7 2 )$ = 47tG(S Pf + 5p Q ) , (10) 

where the perturbed fluid pressure 5pf — WfSpf and 
the perturbations in the quintessence energy density and 
pressure are given by 

S PQ = QSQ-Q 2 ^> + ^SQ, (11) 

Spq = QSQ-QH-^SQ. (12) 

The perturbed Euler-Lagrange equation and fluid con- 
servation equation lead to 
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where Sf = Spt/pt and Vf gives the fluid velocity. From 
now on we use N as a time variable and study the evo- 
lution of <f>, S f , Sq = Spq/pq and Sp = Spq/pq in the 
long- wavelength limit (k/aH « 1). We define the vector 
y = ($, Sf, Sq, Sp) T , and considerable algebra leads to 
the expression 



y' = .F(7f,x) x y , 



(15) 



where the matrix J 7 (7f,x) is given by 
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In order to display the above expression in compact form 
we used the variables Of = 1 — Oq and wq = 7 q — 1, but 
the matrix depends only on the five independent param- 
eters 7 f, Og, 7 q, and 9q. Recall that 7 f and 7 q take 
on values between and 2. 

The general evolution is of course very complicated, 
but as seen in Fig. 1, during each of the five regimes 
described above the coefficients within the matrix take 
on constant values, and this allows us to study the main 
features of the perturbation evolution. We are interested 
in the eigenvalues rii of the matrix T and their corre- 
sponding eigenvectors y^. The solution then takes the 
form 



where the matrix £( 7 f,OQ, 7 g) is given by 



y = ^2 A y, cxp(n 4 A) 
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where the Ai are constants given by the initial conditions. 



A. The adiabatic case 

Before considering the general case, we restrict our- 
selves to adiabatic perturbations. For perturbations to 
be adiabatic, they must share a common perturbation 
according to the prescription 
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which ensures that all matter perturbations vanish 
on uniform-density hypcrsurfaces. Note that the 
quintessence pressure perturbation, as well as its den- 
sity perturbation, must satisfy the adiabatic condition 
(for the perfect fluid adiabaticity of its pressure perturba- 
tion is automatically guaranteed by its equation of state). 
These conditions can be rewritten as 
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It is well known that initially adiabatic perturbations 
remain purely adiabatic, and indeed it is not difficult to 
check that these conditions are conserved through evolu- 
tion by our equations. We can therefore reduce the dy- 
namical system in the adiabatic case to a system of two 
first-order equations. We define the vector z = ($,<5f) T , 
and using Eqs. (15) and (19) we find 
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where 7to t = 7 fOf + 7 qOq- The eigenvalues and eigenvec- 
tors are given in the upper part of Table II. We see, as is 
well known, that there are always a constant and a decay- 
ing adiabatic mode, the former giving the late-time solu- 
tion $ = - 7tot <5 f /2 7f = -5 tot /2, where <5 tot = 6p tot /p tot . 



B. The general case 

We now return to the full set of perturbation equations 
Eq. (15), continuing to consider the regimes in each of 
which the coefficients of the matrix T remain constant. 
We summarize our main results in the lower part of Ta- 
ble II. For the perfect tracker regime the eigenvectors 
y 3 and y 4 have not been given, since they are long and 
complicated formulas which are anyway not very rele- 
vant. In order to simplify some expressions we have used 
the variables 
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z' = £( 7 f,0 Q , 7Q ) x z, 
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For each regime the adiabatic modes can easily be iden- 
tified as the first two entries in the table. Now let us an- 
alyze the two other, non-adiabatic, modes. In the kinetic 
case there are a growing mode, for which Sp = —5q and 
thus SQ — (since $ = 0), and another constant mode 
corresponding to SQ = 0. During the transition and po- 
tential regimes the former growing mode becomes con- 
stant and for each regime the fourth mode is decaying. 
Therefore, before entering the tracker the quintessence 
field may feature large non-adiabatic perturbations. As 
long as the Universe is dominated by the fluid, they are 
isocurvature perturbations. 

In the usual tracker case the last two eigenvalues may 
have an imaginary part, leading to oscillations, and their 
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FIG. 2: Features of the eigenvalues nz and n 4 in the case of 
an usual tracker regime, as a function of 7f and 7q. Their 
real part is negative below the solid line and they have an 
imaginary part on the right-hand side of the dashed line. 



more slowly than pf, and hence jq < jf. As one can 
see in Fig. 2, this implies Re(n 3i 4) < . In the perfect 
tracker case one easily sees that the last two modes de- 
cay, possibly oscillating as they do. As long as all the 
other modes are decaying, during the tracker regime the 
constant adiabatic mode y 4 is an attractor. As a re- 
sult, a long tracker period implies the suppression of all 
non-adiabatic modes [5, 6]. Moreover, in the case of a 
sub-dominant tracker the late-time evolution of the per- 
turbations is even independent of the quintessence field 
initial conditions. 

Finally, during the domination regime, which is 
reached after the present epoch, the non-adiabatic modes 
are constant and decaying. 



IV. CONCLUSIONS 
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TABLE II: Top: eigenvalues and eigenvectors of the matrix 
Q (adiabatic case). Bottom: eigenvalues and eigenvectors of 
the matrix T according to the different regimes; n ut i, n ut2 , 
n u t3, TCpti and n pt2 are given in Eqs. (22)-(26). 



real part can either be negative or positive according to 
the value of 7f and jq. The regimes of the eigenvalues 
are shown in Fig. 2. However, since the quintessence field 
has to dominate at the present epoch, pQ must decrease 



We have derived four first-order equations to describe 
large-scale perturbations in quintessence scenarios. Dur- 
ing each of the five main regimes of quintessence be- 
haviour, these equations have constant coefficients, en- 
abling analytic solution of the perturbations by eigen- 
vector decomposition. 

We have seen that during the kinetic regime there is a 
growing isocurvature mode which then remains constant 
until tracking begins. However, if the quintessence field 
undergoes a long period of tracking, there remain only 
adiabatic perturbations which, in case of a sub-dominant 
tracker, are independent of its initial conditions. A low 
initial quintessence energy density [8], or a long kinetic 
period, may prevent the non-adiabatic modes disappear- 
ing completely. 

It is possible in principle to carry out the same analy- 
sis without using the large-scale approximation, although 
the equations may then be too complicated to be useful. 
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